1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
| library(MASS) library(class)
x0<- MASS::mvrnorm(n=10, mu = c(1,0), Sigma = diag(2)) x1<- MASS::mvrnorm(n=10, mu = c(0,1), Sigma = diag(2))
make.data<-function(x0,x1,n,m,Sigma1 = diag(2)/5,Sigma2 = diag(2)/5){ r<-length(x0[,1]) x<-matrix(0,n+m,2) for(i in 1:n){ s<-sample(1:r,1) x[i,]=MASS::mvrnorm(n=1, mu = x0[s,], Sigma = Sigma1) } r<-length(x1[,1]) for(i in 1:m){ s<-sample(1:r,1) x[i+n,]=MASS::mvrnorm(n=1, mu = x1[s,], Sigma = Sigma2) } y<-matrix(c(rep(0,n),rep(1,m)),n+m,1) data<-data.frame(y,x) return(data) } train.data<-make.data(x0,x1,100,100,diag(2)/5,diag(2)/5)
|